function [output] = apply_lognormal (variance, input)
% Function for lognormal channel effects

sigma = sqrt(variance);
len = length(input);
output = rand(1,len);
ran = randn(1, len).*sigma;

for i=1:len
    output(i) = exp(ran(i))*real(input(i));
end

 